home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / PROGMISC / PCSSP.LZH / PC-SSP.ZIP / POLYSPTP.ZIP / CNP.FOR next >
Text File  |  1985-11-29  |  2KB  |  55 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE CNP
  5. C
  6. C        PURPOSE
  7. C           COMPUTE THE VALUES OF THE CHEBYSHEV POLYNOMIALS T(N,X)
  8. C           FOR ARGUMENT VALUE X AND ORDERS 0 UP TO N.
  9. C
  10. C        USAGE
  11. C           CALL CNP(Y,X,N)
  12. C
  13. C        DESCRIPTION OF PARAMETERS
  14. C           Y     - RESULT VECTOR OF DIMENSION N+1 CONTAINING THE VALUES
  15. C                   OF CHEBYSHEV POLYNOMIALS OF ORDER 0 UP TO N
  16. C                   FOR GIVEN ARGUMENT X.
  17. C           Y     - RESULT VALUE
  18. C                   VALUES ARE ORDERED FROM LOW TO HIGH ORDER
  19. C           X     - ARGUMENT OF CHEBYSHEV POLYNOMIAL
  20. C           N     - ORDER OF CHEBYSHEV POLYNOMIAL
  21. C
  22. C        REMARKS
  23. C           N LESS THAN 0 IS TREATED AS IF N WERE 0
  24. C
  25. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  26. C           NONE
  27. C
  28. C        METHOD
  29. C           EVALUATION IS BASED ON THE RECURRENCE EQUATION FOR
  30. C           CHEBYSHEV POLYNOMIALS T(N,X)
  31. C           T(N+1,X)=2*X*T(N,X)-T(N-1,X),
  32. C           WHERE THE FIRST TERM IN BRACKETS IS THE ORDER,
  33. C           THE SECOND IS THE ARGUMENT.
  34. C           STARTING VALUES ARE T(0,X)=1, T(1,X)=X.
  35. C
  36. C     ..................................................................
  37. C
  38.       SUBROUTINE CNP(Y,X,N)
  39. C
  40.       DIMENSION Y(1)
  41.       Y(1)=1.
  42.       IF(N)1,1,2
  43.     1 RETURN
  44. C
  45.     2 Y(2)=X
  46.       IF(N-1)1,1,3
  47. C
  48. C        INITIALIZATION
  49.     3 F=X+X
  50. C
  51.       DO 4 I=2,N
  52.     4 Y(I+1)=F*Y(I)-Y(I-1)
  53.       RETURN
  54.       END
  55.